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Several recent experiments have reported an anomalous temperatnre dependence of the Coulomb 
drag effect in electron-hole bilayers. Motivated by these puzzling data, we study theoretically a 
low-density electron-hole bilayer, where electrons and holes avoid quantum degeneracy by forming 
excitons. We describe the ionization-recombination crossover between the electron-hole plasma and 
exciton gas and calculate both the intralayer and drag resistivity as a function of temperature. The 
latter exhibits a minimum followed by a sharp upturn at low temperatures in a qualitative agreement 
with the experimental observations [see, e.g., J. A. Seamons et al., Phys. Rev. Lett. 102 , 026804 
(2009)]. Importantly, the drag resistivity in the proposed scenario is found to be rather insensitive 
to a mismatch in electron and hole concentrations in sharp contrast to the scenario of electron-hole 
Cooper pairing. 

PACS numbers: 71.35.Ee, 73.63.Hs 


Coulomb drag effect is a sensitive probe of interactions 
and collective phases in bilayer systems (see Ref. Eli 
for a review). In its usual setup, an electric current in the 
first layer, /drive; drags charge carriers in the other one. If 
the second layer is closed, the drag force is compensated 
by the Coulomb force induced by a voltage drop, Vdrag; 
and the drag resistivity of the bilayer pD = bdrag/Zdrive 
is measured. If the bilayer involves two weakly-coupled 
Fermi liquids, the temperature dependence of the drag re¬ 
sistivity at low temperatures is quadratic pu T^, which 
is well established both theoretically and experi¬ 

mentally [i,0- Any deviations from that Fermi-liquid 
behavior can signal the appearance of collective phases 
or correlations in the bilayer system. 


A number of recent experiments (S-Ul on the electron- 
hole CaAs/CaAlAs bilayers have observed an anoma¬ 
lous temperature dependence of drag resistivity at at 
the intermediate doping n^^h) ~ 5 10^° cm“^. The T- 
dependence of pD was shown to achieve a minimum, fol¬ 
lowed by a growth and saturation at lower temperatures, 
which were rather insensitive to the concentrations mis¬ 
match (see also related experiments for other realizations 
of electron-hole bilayers EMI). This behavior cannot 
be attributed to interlayer exchange and correlation ef¬ 
fects 15l-ll7j|. that are relevant in that regime, and does 
not appear for electron-electron and hole-hole bilayers for 
similar parameters. Therefore, there is strong evidence 
for an excitonic origin of the effect, but its detailed un¬ 
derstanding is still lacking. 


There were a number of theoretical attempts to explain 
the experiments based on the Bardeen-Cooper-Schrieffer 
(BCS) model of electron-hole Cooper pairing 18l-21|. 
which is valid in the high doping re gim e and can be 
the origin of the dipolar superfluidity |18l . [il . . The 
mean-field theory predicts a jump of drag resistivity at 
the pairing temperature to a value comparable to a sin¬ 
gle layer resistivity [1^. The jump can be smoothed by 
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FIG. 1: (Color online) Shown is the dependence of the 
drag resistivity pD on temperature for matched concentra¬ 
tions of electrons and holes with excitons [Fig. 1(a)] and 
without excitons [inset 1(b)]. The curves correspond to 
different density per layer denoted by their values in 

cm“^, and axes of the inset coincide with ones of the main 
plot. The drag resistivity pu achieves a minimum at Tu 
within ionization-recombination crossover between the high- 
temperature regime, T ^ Aexc, where the drag is dominated 
by Coulomb interactions in the electron-hole plasma, and the 
low-temperature regime, T « Aexc, where the drag is domi¬ 
nated by excitons. 


pairing fluctuations [lElii, which are a precursor to the 
paired state, and both Aslamazov-Larkin 27| and Maki- 
Thomson contributions are important here. How¬ 

ever, Cooper pairing and the fluctuations are very sensi¬ 
tive to the mismatch 31, ^ in contrast to experimental 
observations. 


Here we present an alternative theoretical scenario for 
the effect involving the formation of excitons, which are 
a bound state of spatially separated electrons and holes, 
with a small binding energy, /?exc- For T /?exc; exci¬ 
tons ionize to form a classical electron-hole plasma and 
the drag effect is dominated by the Coulomb interactions. 
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At low temperatures, the appearance of excitons strongly 
enhances the drag and single-layer resistivities, leading 
to the upturn in the former. The anomalous behavior 
is robust against the mismatch in the concentration of 
electrons and holes: while the magnitude of the upturn 
is affected by it, the temperature Td, where the drag re¬ 
sistivity reaches minimum is insensitive to the mismatch. 
Proposed scenario is valid and self-consistent at low dop¬ 
ing, and the calculated excitonic upturn is considerable 
larger than the observed one. Nevertheless, our results 
are in a qualitative agreement with the existing exper¬ 
iments. The main conclusion of our work is that the 
picture of exciton formation is more relevant to the inter¬ 
mediate doping regime in experiments, than the scenario 
of electron-hole Cooper pairing. 

Model and the excitonie crossover. The system of spa¬ 
tially separated electrons and holes, which can bind to 
form excitons, can be described by the following Hamil¬ 
tonian 


H = E ^exc (p)^pS^p5' + 

pS 


ps 


4-1 T/ /' 'l + + 


pp q 

ss' OLOi' 


Here Uaps and &ps are annihilation operators for elec¬ 
trons (a = e = 1), holes (a = h = —1) and exci¬ 
tons with momentum p and internal degeneracy spin 
index s = (| i), | t)) and 5 = (| U), I tJ-)), I it), I tt))- 
Their dispersions are eQ,(p) = /2ma and eexc(p) = 
p^/2mexc-^^exc with TTiexc = me-l-mh and Aexc being the 
exciton mass and its binding energy; 14(a(q) = 27re^/eg 
and T4<a(q) = —27re^e“'^'^/eg are bare intralayer and in¬ 
terlayer Coulomb interactions with interlayer spacing d 
and bare dielectric permittivity e. We do not specify 
the interaction with disorder explicitly, but assume re¬ 
laxation times Ta and Texc to be momentum indepen¬ 
dent, which implies the short-range disorder to be the 
dominant scattering mechanism. 

For all numerical calculations we use the set of pa¬ 
rameters related to the GaAs/GaAlAs bilayer in exper¬ 
iments TUe Ri 0.067mo, mh ^ 0.4mo, d ~ 30nm, 
e = 12.4 with uiq the bare electronic mass. The re¬ 
laxation times Ta are parametrized by mobilities Me Ri 
2 10® cm^/Vs, and Mh R! 3 10® sm^/Vs. The exci¬ 
tonic relaxation time, Texc = TO*re7ii/(reme -I- Thmh), 
where the reduced mass is m* = meTOh/(me -I- mh), 
corresponds to the mobility Mexc R^ 3.4 10^ cm^/Vs. 
Nevertheless, excitons being nonlocal objects are more 
sensitive to interlayer tunneling and other factors, so 
their mobility can be considerably reduced and here 
we use Mexc ~ 10^ cm^/Vs. The effective Bohr ra¬ 
dius, Ob = r; 11.8 nm, and Rydberg energy, 

Eb = m*e^/2fi^e^ = 55.4 K, give the spatial and energy 
scales. The exciton energy, Aexc, can be considerably 



FIG. 2: (Color online) The temperature dependence of con¬ 
centrations of electrons Ue, holes rih and excitons riexc, which 
are given by the Eqs. (H), for fixed total concentrations per 
layer n® = 8 10®cm“^, — 6 10®cm“^. The dependen¬ 

cies for other values of Ue(h) arc qualitatively similar. At 
high temperatures T Aexc there is a long excitonic tail 
riexc ~ Werih/n* ~ T~^. The concentration of excitons at zero 
tempearture is equal to that of the minority species (holes in 
this case) in the limit of large temperatures. 


smaller than Eb at d > ob and is sensitive to screening, 
so here we use Eexc ~ 0.5 K, corresponding to the ex¬ 
citon size Oexc Ri 110 nm, as an independent parameter. 
The model is self-consistent if excitons weakly overlap, 
which corresponds to the doping ne(h) ^ 10^° cm“^. 

The ground state of the model is believed to be the 
exciton condensate that forms at the temperature Tq < 
Eexc and can coexist with the degenerate gas of electrons 
or holes in the presence of their concentration mismatch. 
However, below we focus on the ionization-recombination 
crossover regime T > Eexc, where the distributions of 
electrons, holes and excitons are non-degenerate. To 
calculate their concentrations we recall that in experi¬ 
ments the total concentrations of charged particles per 
layer n® are controlled independently by electrical dop¬ 
ing, so riexc + na = ria- Here Ua and Uexc are concentra¬ 
tions of quasiparticles. Reintroducing the grand canon¬ 
ical Hamiltonian, EIq = El — J2a Eainexc + fia), with 
chemical potentials pa as Lagrange multipliers, we get 
the chemical potential of excitons as pexc = Me + Mh- 
The equation for concentrations can be simplified to 
nenh/n-» + Ua = n®, where the concentration n* = 
m*T exp[—Aexc/r]/(27rM^). The temperature dependen¬ 
cies of fermionic and excitonic concentrations are given 
by 


^exc — 


Sn° - n* 


-I- u* 


((5nO)2 + nl+ 2n*nf), ; 
((5n0 )2 + nl + 2n*n!^ , 


( 2 ) 


where Sn^ = n® — n® and are the concern 
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tration mismatch and the total concentration. 

The temperature dependence of the concentrations is 
depicted in Fig. 2. At low temperatures T Aexc the 
fraction of unbound electrons and holes is exponentially 
small, while within the crossover T > E^xc there is a 
long non-degenerate tail of excitons decreasing as T~^ 
according to riexc ~ rign^/n^:. The exciton gas can be 
considered non-degenerate until Tq rs 0.3 K [33| . 

Phenomenology of the drag effect. In the presence of 
electrons, holes and excitons the conductivity tensor of 
the bilayer system is given by 

f Je \ _ f ^exc “t” f^exc f Eq 

\ Jh J \ ^exc f^exc E J \ Ah 



where ctq, = nae^Ta/ma and cTexc = nexce^Texc/w-exc are 
their Drude conductivities. Excitons, being composed of 
electrons and holes from different layers, contribute to 
both diagonal and off-diagonal components of the con¬ 
ductivity tensor with opposite signs. The transconduc¬ 
tivity ctd originates from the Coulomb interaction be¬ 
tween electrons and holes and is calculated microscopi¬ 
cally below. The drag resistivity pD and single layer re¬ 
sistivities Pa, being the components of the inverted con¬ 
ductivity matrix ([3]), can be written in a compact way 


PD{a) — 


0'D(o 


(4) 


(TeCTh + {(Je + 0-h)cr exc 

At zero temperature the excitonic contribution domi¬ 
nates and they become 


PD 


(n° - n%)e^Ta 


Pa — Pd + 


©QQ,?7T,exc 

^“e^rexc 


(5) 


Here Qaa = E){na — ria) is the Heaviside function. If 
densities of electrons and holes are perfectly matched, 
both single layer resistivities Pa and pd diverge at T = 
0. This corresponds to an insulating excitonic ground 
state with the perfect drag effect: the relation between 
the electric current in a layer, induced by a current in 
the other layer, is /drag = —/drive- Our considerations 
assume T 3> Tq, where there is a competition between 
CTexc and CTD, but the zero-temperature values ([5]) reflect 
the strength of the low-temperature upturn. 

Electron-hole transconductivity. The transconductiv¬ 
ity ctd can be calculated in the second order of perturba¬ 
tion theory in the interlayer Coulomb interaction as 
follows 


1 


Ctd = 


IGttT 


E 


/ . ,2, u> J xe (q,W,w)x 

-oo Smh {^) (0) 

riI'(q,Cd,cd)|C/eh(q,w)^ 


where U{q,uj) is the screened interlayer interaction and 
r^^(q, WjW) is the current-charge-charge nonlinear sus¬ 
ceptibility. If the relaxation times Ta are momentum in¬ 
dependent, as we assume here, it is given by 

r^a n^2(q:‘^), (7) 

rUa 


where n^ 2 (qi is the imaginary part of the polarization 
operator, which for a non-degenerate gas is given by 


H 


R _ 

a2 — 


Tq 


sinh 


r uj 1 


r 2 21 

g> g 

[2t\ 

exp 

4T2g2 4/2 _ 


( 8 ) 


Here the fa = is the characteristic thermal mo¬ 

mentum scale. For the interaction [/(q, w), the static 
Debye-Hiickel approximation, that ignores the presence 
of neutral excitons, yields 


Aeh(q) = 


27re"^ 


qe 


— qd 


{q + Ke){q + Kh) - KeUhe 


(9) 


Here «:a(q) = K^f,^{q/qa), = 2Tre‘^na/eT is the 

Debye-Hiickel screening momentum and f{x) is the di¬ 
mensionless function /(2a;) = v^exp[—a;^]Erh(a;)/2x 
with Erfi(a:) to be the imaginary error function. The 
static screening approximation does not take into account 
possible plasmon contribution 34, [s^ . which consider¬ 
ably enhances the drag effect for 0.4 ^ T/p < I. How¬ 
ever, in the non-degenerate regime, the plasmons become 
strongly damped, and can be ignored. The integral over 
frequencies in Eq. ([5]) can be calculated explicitly and we 
get 


/— 2 4 

^/^^e TeTh gd r 

Ctd = —tt:- Jo 


32 h h? rUf 


rrih 


( 10 ) 


with momenta qa = hd g* = \/2m*T and a dimen¬ 
sionless integral Iq over rescaled momentum q given by 


4 = 


dx 




0.0^3 


C — 4q^ d'^ 


[{x + Ked){x + Klid) — Ke Kiid'^e 


( 11 ) 


There are three different momenta gd, g*, (for calcu¬ 
lations of asymptotes we assume that and have 
the same order of magnitude) in the integral Iq, and 
the characteristic momentum, transfered between elec¬ 
tron and hole layers, is the smallest of them. If these 
momenta are well separated, the asymptotic behavior of 
the integral Iq can be evaluated analytically. There are 
four different regimes: / : g* <C gd, //-i- : gd -C g*, ; 

//_ : < gd and III : < gd < g* with 


/: I, 




11+ , Iq 


//_ : Iq = y/f^qlnlnldT, III : Iq 


q^,d ^ 

g»«:g/tj(ci^ 

2 ■ 


Regimes I {T < T^) and III {T^ < T) appear at 
small and large temperatures. Depending on the con¬ 
centration ne(h) one of //_ and 11+ is between them 
{T^ < T < T^). The corresponding boundaries are 
given by T 2 = iTrEBnaOBd, Tff = Tj^ = /?B(aB/d)^ and 
Tf = EBi^nnaa^)"^^^. The point at which T)^ = 
and the regimes II± merge and disappear corresponds to 
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FIG. 3: (Color online) (a) and (c) The temperature dependen¬ 
cies of the drag resistivity pD and the resistivity of electrons 
Pe in the presence of the mismatch in electron and hole con¬ 
centrations. (b) and (d) The corresponding values at zero 
temperature, which are given by Eqs. ©)• The strength of 
the excitonic enhancement of both pD and p^ is defined by 
the mismatch, while the temperature dependencies are quite 
insensitive to it. 


ni 2 = QB/'i'n'cP = 3.1 10® cm“® and T 12 = EB{aB/d)'^ ~ 
7.1 K. For the densities of interest, the momentum scales 
are not well separated, the range of the applicability of 
the asymptotes (fT^ is reduced to T <C and T T^. 
Below, we calculate Iq numerically. Drag resistivity of the 
bilayer. First, it is instructive to analyze the dependence 
of drag resistivity pD on the temperature T ignoring the 
presence of excitons. It is shown in Fig. 1-b (for matched 
concentrations of electrons and holes). At high tempera¬ 
tures T 2 ^ T, the screening disappears, making k® the 
smallest momentum scale, and the drag resistivity de¬ 
creases as pd jd. In the intermediate regime 

^ T , the scattering momentum is and 

the asymptotic form is pD In'ln^df'. These two 

scattering regimes are usual for bilayer fermion systems 
(along with regimes where plasmons [341 l36l| and phonons 
[3711 dominate and the hydrodynamic one [381139|), but 
the latter corresponds to pD ~ T® due to degeneracy of 
fermions. For the considered system, at low tempera¬ 
tures T < , the electrons and holes avoid degeneracy 

by transforming into excitons and their characteristic mo¬ 
mentum scale g* becomes the scattering one leading to 


the asymptotic behavior pD ~ That regime 

usually does not appear in a fermionic bilayer due to 
quantum degeneracy of fermions. 

The temperature dependence of pD in the presence of 
excitons are presented in Fig. I-a (perfectly matched den¬ 
sities) and Fig. 3-a (with a mismatch). The latter is sup¬ 
plemented by the inset Fig. 3-b in which the dependence 
of Pd on the mismatch at zero temperature is depicted. 
The long excitonic tail, which weakly depends on tem¬ 
perature, considerably enhances the drag resistivity even 
at high temperatures T ^ Aexc- The dependence has a 
clear minimum at the temperature Td, which lies within 
the crossover Flexc ^ ^ T^- The strength of the up¬ 

turn is defined by the mismatch, while the temperature 
TB,{n°,n^) smoothly increases with both its arguments 
and does not have any features for the matched case. 
This makes the minimum in the temperature dependence 
of Pd shallower with increasing of both concentrations, 
as seen in the experiment. The excitonic contribution to 
the drag resistivity pD can be well-fitted by a combina¬ 
tion of functions T~^ and T“®. The former dominates 
at high temperatures T ^ Flexc, while the latter plays 
the major role at T ~ Aexc- At lower temperatures the 
drag resistivity saturates to a value, which depends on 
the imbalance of concentrations (see Fig.3-b). 

Resistivity of electrons is presented in Fig.3-c and sup¬ 
plemented by the inset (d), where its dependence on the 
mismatch at T = 0 is depicted. Depending on the mis¬ 
match, its enhancement vary by an order of magnitude, 
while the temperature dependence is quite insensitive to 
it. The dependencies for the resistivity of holes are qual¬ 
itatively the same. 

Discussion. The proposed scenario of genuine exci¬ 
tonic drag effect does not assume any phase transition 
and/or coherence of excitons, which in our model may 
occur at lower temperatures Tq (Localization effects and 
their interplay with other ground states, not involving 
exciton condensation, can not be ruled out: e.g., an 
excitonic Bose glass [i^ |4lj or an exotic Bose-metal 
phase [ 4 ^, which was conjectured to exist in models 
involving dirty composite bosons and gapless fermionic 
excitations). We argue however that the upturn in pD 
is unrelated to the quantum effects including localiza¬ 
tion, but appears at the temperature Td correspond¬ 
ing to the ionization-recombination excitonic crossover 
Aexc ^ 2b ^ from a classical electron-hole plasma 
to a classical exciton gas. The exact value of the Td 
is non-universal and depends on the interlayer distance, 
quasiparticle mobilities, effective masses, etc. 

For explicit calculations above, we have used a range of 
electron and hole concentrations, which is about an order 
of magnitude smaller than the ones in the published ex¬ 
periments to ensure that the assumptions of our model 
are self-consistent. In the intermediate doping regime 
realized in experiment so far, the excitons overlap and 
cannot be considered as two-particle objects anymore. 
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To develop a quantitative many-body theory for the 
Coulomb drag effect in this intermediate regime is diffi¬ 
cult, because of complicated interplay of the Pauli block¬ 
ing effects, self-consistent screening, and coexistence of 
excitons with degenerate gas of electrons and holes. The 
extrapolation of our results to this regime considerably 
overestimates the strength of the excitonic upturn seen 
in experiments. Nevertheless, the observed behavior of 
the drag resistivity on temperature and concentrations is 
qualitatively captured, and we conclude that the picture 
of exciton formation is more relevant to the experiments, 
than the scenario of electron-hole Cooper pairing and 
pairing fluctuations. 
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Supplemental Material: ’’Anomalous Coulomb Drag in Electron-Hole Bilayers due to 
the Formation of Excitons” by Dmitry K. Efimkin and Victor Galitski 


The Supplemental Material presents a calculation of the polarization operator n(q, w) of two-dimensional nonde¬ 
generate electron gas. The general expression for the polarization operator n(q, w) is given by 


n(q,w) = 2^ 

P 


^f(Cp) 

-I- a; 


^F(^p+q) 

■Cp-l-q T iS 


(SI) 


Here = p^/2m — /i is the dispersion law of fermions and nF(fp) ~ exp[—/3^p] is their Maxwell distribution function 
with chemical fi = T\n(Trh^n/mT), where n is the total concentration of fermions including spin degeneracy. It is 
instructive to calculate the imaginary part of the polarization operator n 2 (q, w) at first and it can be presented as 
follows 


U 2 {q,w) = -TT—(1 - e" 


0^ / 


pdp _ 
e ' 
m 


r2TT 


d(j) ( r pq 

- TT- - — cos {4>) 


27T 


2m 


m 


(S2) 


The character energy and momentum of nondegenerate electron gas are wt = T and gx = V^mT, and the polarization 

operator, which has the dimension of the density of states, has the scale n/T. As a result, introducing dimensionless 

units p = qT\/x, q = qrq, w = wtw, and n(g, uj) = —nIL{q, <D)/T we get 

n 2 (g, w) = 7r(l — e“‘^) f dxe~^ f (w — — 2g-v/x cos ^) . (S3) 

Jo Jo 27r 

The argument of the delta function achieves zero only if x > xo, where xq = (w — g^)^/4g^, that leads to 


n2(g,w) = 


4g 


(ix^2_ J j ^ _ cos(^) ) = 


(l-e~^) 

2g 


dxe 


\/x - xo 




2g 


-, (S4) 


and finally we get the expression ([5]) from the main part of the paper, written in dimensionless units 


n 2 (g,a;) = — sinh 

q 


/UJ\ 

■ g2- 

(2) ""P 

_ 4g2 4_ 


(S5) 


The dependence 112(g,d)) on momentum and frequency is presented in Fig. Sl-a. The value 142(g,w) is nonzero for 
arbitrary momentum and frequency and has prominent maximum at tc ~ 1 and g ^ 1, which correspond to the 
thermal frequency ujt and momenta gx. The calculated dependence for nondegenerate electron gas drastically differs 
from one at zero temperature. The value n 2 (g, w) at zero temperatures is nonzero only within the continuum of 
two-particle excitations, which has sharp boundaries. 

In the paper we use the static screening approximation and need only real part of the polarization operator IIi(g, 0) 
in the static limit. It can be calculated from the imaginary part with help of Kramers-Kronig relations as follows 


ni(g,0) 


duj n2(g,a;) 
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3g2 — 2qi + 6 


(S6) 


Here Erfi(g/2) is the imaginary error function. The approximate expression interpolates the corresponding asymptotics 
and well describes the momentum dependence of ni(g, 0) at arbitrary momentum. The momentum dependence of 
ni(g, 0) is presented in Fig. Sl-b. 
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FIG. SI: a) Frequency and momentum dependence of 
dependence of the real part of the polarization operator 



le imaginary part of polarization operator 112 (g,a)); b) Momentum 
i{q,Ld) in the static limit. 









